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Abstract 

We consider a particular 4 state spin system composed of two Ising spins ( Sj,, ax ) 
with independent hopping parameters ki,K2, coupled by a bilinear Yukawa term, 
ysxax- The Yukawa term is solely responsible for breaking the global Z2 x Z2 sym¬ 
metry down to Z2. This model is intended as an illustration of general coupled Higgs 
system where scalars can arise both as composite and elementary excitations. For the 
Ising example in 2 d, we give convincing numerical evidence of the universality of the 
two spin system with the one spin Ising model, by Monte Carlo simulations and finite 
scaling analysis . We also show that as we approach the phase transition, universality 
arises by a separation of low mass spin waves from an extra set of spin waves with an 
energy gap that diverges as the correlations length diverges. 


1 Introduction 


It often happens in quantum field theory that exactly the same physical system can be 
described by very different actions. For example, the critical behavior of a two dimensional 
paramagnetic him can be described either by a Bosonic 2d Ising model, or equivalently 
by a 2d free Fermion held []^. Another example is the exact equivalence of Sine-Gordan 
theory and the massive Thirring model. Kinks in the elementary scalar held theory become 
elementary Fermions while the scalars are composite Fermion anti-Fermion bound states[§]. 

A similar equivalence exists for the Goldstone modes of chiral symmetry breaking which 
can be described either in terms of an elementary scalar held theory or as composite scalar 
modes in a Fermion-antiFermion channel. For example it has been shown that the Higgs- 
Yukawa model and the generalized Nambu-Jona-Lasinio Model are identical. Note that 
in the last example, the number of degrees of freedom appears to be very diherent at the 
cut-oh scale, yet the low energy physics turns out to be exactly the same. All of these are 
manifestations of the the well known universality property: diherent microscopic Lagrangians 
can lead to the same continuum theory (critical behavior) as long as the hxed point structure 
and the low mass excitations of the two actions are the same. 

Here we wish to consider a diherent but closely related phenomenon where the same 
excitation is represented twice in a single Lagrangian by two diherent helds, yet in the 
scaling limit the theory is universal to either spin system taken separately. Such a redundant 
description has recently been advocated as an improved lattice formulation of GGD in the 
so called chirally extended action (XQGD)[0,p. Here an elementary scalar held is coupled 
to the Wilson lattice QGD action. It was conjectured that XQGD is in the same universality 
class as the Wilson lattice QGD action and its critical property describes the continuum 
QGD physics. An obvious problem for XQGD is the extra degrees of freedom introduced by 
the scalar helds on the lattice. Will the extra degrees of freedom decouple in the continuum 
limit? This question can not be answered easily without nonperturbative calculations. 

An important feature of this suggestion is the existence of a global (Higgs) symmetry G 
in each held theory so that the combined theory has & G x G symmetry that is broken down 
to G by a Yukawa coupling. In these more complicated examples, it is not obvious when they 
remain in the same universality class as the original theory or if numerical methods can give 
convincing evidence. We were able to prove that an extended Nambu-Jona Lasino model 
and a two spin coupled 0(N) models are equivalent in large N Si- At large N, the extra 
degrees of freedom become decoupled in the continuum limit and the low energy spectrum is 
exactly the same as if no extra degrees of freedom were introduced at the cut-oh scale. Our 
purpose in this article is to begin the use of numerical methods to supplement analytical 
methods. To this end we consider a simple model of redundant coupled Ising spins to see if 
Monte Garlo methods can settle the question of universality and expose the mechanism by 
which it is achieved in the scaling limit. 
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2 Yukawa Coupled Spin Ising Model 


For the our two spin Ising system, we consider the action, 

S ^ ^ ^2 ^ ^ ^x-\-fi^x y ^ ^ ^x^x 1 ( 1 ) 

X^[l X,fl X 

where = ±1 and = ±1- The sum over x = {ni,n 2 ) is performed on a two dimensional 
lattice, rij = 1, 2, • • •, L, with nearest neighbor couplings, fi = (1, 0), (0,1). A very general 
class of four spins models could be considered. Our particular choice we call a “Yukawa” 
coupling borrowing terminology from held theories for spontaneous breaking of continuous 
symmetry. When the Yukawa coupling y is set to zero, this model reduces to two 
independent 2d Ising models with a global Z 2 x Z 2 symmetry. But when y is different from 
zero, only a single Z 2 symmetry is survives. This is a special example of 4 state spin models, 
which to our knowledge has no analytic solution. It is no longer self-dual. 

At the cut-off scale there are two helds, ax and Sx, in the coupled Ising model in Eq. (|l|). 
However, as we will show in the following sections, only one eigenvalue of the mass matrix 
remains in the spectrum in the continuum limit. The other mass eigenvalue stays at the 
cut-off scale. We believe that the light mass corresponds to the massive particle of a 2d Ising 
model. We will further show that the coupled Ising model follows the scaling laws of a 2d 
Ising model. Therefore, our model in Eq. (|l]) is in the same universality class of a 2d Ising 
model. 

Admittedly in this case the universality of our 4 state spin model to the Ising model is 
not surprising. The Yukawa coupled two spin model is in fact nothing but an Ising model 
written on 3-d lattice with only two layers in the third axis and three independent couplings: 
Ki and K 2 in each x-y layer and y along the z axis. Since this model can be viewed as a 
d = 2-|-e model with vanishing e in the thermodynamic limit and since other anisotropic Ising 
models are known to be in the isotropic Ising universality class, universality with a single 
Ising spin system appears likely . When we tune to the critical surface and the correlation 
length S, becomes much larger than 2, the system should behave like a single layer 2d Ising 
model. Still it is important for our larger purposes to see if numerical methods can give 
convincing evidence and that analysis of the spectrum can illuminate the mechanism for 
approaching the single spin system in the scaling limit. Finally we should note that since 
the totally decoupled system {y = 0) is clearly not a single Ising system, so one expects that 
by carefully tuning the Yukawa term a new interacting two spin fixed point must also exist 
which is not universal to a single spin system. This fine tuning limit is what is often required 
to obtain a true 2 Higgs system in the continuum. 


3 Phase Diagram 


First we need to determine the phase diagram of the model. For arbitrary ki, k ,2 and ?/, Eq. 
(|^ describes a two layered Ising model with different couplings in each layer. In general. 
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we expect that the magnetization in one layer will induce the magnetization in the second 
layer when the coupling between layers, ?/, is nonzero. Therefore, the size of the symmetric 
region,as shown in Figure |^, will be reduced. A few limiting cases can be easily derived. 

If ?/ 7^ 0, ^2 = 0, Hi 7 ^ 0, the sum over can be carried out explicitly 


Z 


EEe 


X/Xj/i ^x^x+fi T y X/a; 

{A 

x(^x+ii^^ ln[ch{ya^)] 

const X ^ '^x,y. C^xC^x+y ^ 

{o-} 


( 2 ) 


(3) 


Thus the model becomes the standard 2d Ising model and the critical point is known from the 
exact solution Hi^cr = Hc = 0.44068.... It is also easy to show that a spontaneous symmetry 
breaking in a will induce a magnetization in s 

< s >= tanh{y) < a > , (4) 

where < a > is given by the 2d Ising model solution. Similar result can be derived for 

Ki = 0, K2 7^ 0,2/ 7^ 0. 

In the limit y ^ oo, would be completely aligned with Sx, and up to a constant the 
action can be written as 

S = -{Ki + K2)'^(Jx(Jx+f, . (5) 

This is again a 2d Ising model with critical point given by a straight line. 


Ki + K2 = He = 0.44068... , 


( 6 ) 


in the Hi — H 2 phase plane. At hnite y, the phase boundary lies between the limit of ?/ = 0 
and ?/ —> oo as shown in Figure ^ It was determined using Monte Carlo simulations which 
were carried out on a 64-node CM-5 at Boston university. For simplicity, we have used 
the heat-bath algorithm. The size of lattices we used range from 16^ to 128^. The phase 
boundary was determined by observing the magnetization, dehned as: 


< mi >=< 



>, 


< m2>=< 



> 


(7) 


It is non-vanishing on a hnite lattice in both phases, but will converge to the physical 
magnetization in the inhnite volume limit. Near the phase transition point, < mi >, < m 2 > 
changes rapidly from hnite values to almost zero in the symmetric phase. This cross-over 
behavior becomes sharper with increasing lattice size. The results for y = 0.50 are shown 
in Figure The points have been computed while the line is a simple interpolation. /,From 
this hgure one can get that y = 0.50, H 2 = 0.20 and Hi = 0.325(5) is a critical point. This 
point will be used in the study of the propagators. For a more precise measurement of the 
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Figure 1: Phase diagram for y = 0.50. The lattice size is 32^. The solid line is an 
interpolation to the points. 

critical point, we used the Binder index |P] 


Ui{L) 


< mf > 
3 < mf 


U2{L) 


< > 
3 < ml 


( 8 ) 


For the simplest case ki = ^2 = 2//2 = /9, we have Ui = U 2 = U{L, k). When the volume 
is large, U{L, n) should cross at an unique point (f/*, Kc). However, for practical L and L', 
the crossing point depends on the ratio s = L'f L due to residual hnite size corrections. We 
obtained the crossing point for a series of values of s, and htted to the form ^ 

- = — _^ — (9) 

P ln{s) Pc ’ 

as shown in Figure We found Pc = 0.2759(2). Another way to locate the critical point is 
to use the hnite size scaling behavior of the susceptibility 0 |^. We htted the location of 
the peak (Pmax) of susceptibility to: 


Pv 


1 

Wc 


( 10 ) 


with iz = 1 (see next section). We found (see Figure |) Pc = 0.2766(4), which is consistent 
with the result obtained using the Binder index. To be conservative, we take the diherence 
in the above two estimates as an indication of the true error and quote the critical point to 
be pc = 0.2759(7). 
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Figure 2: Determination of the critical point using the finite size scaling of the Binder 
Index. The fitting curve is l/fi = a/s + 1/fic (s = log{Ll/L2)). The computed values of the 
fitting parameters are: fie = 0.2759(2) and a = 0.000(2). 



Figure 3: Location (fimax) of the Susceptibility Maximum vs. L Using Eg\T^ 

(I / fimax 

1/ fic + ajL~'' ) we compute the critical point fic = 0.2766(4). (o is assumed to be u = 1) 
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As an aside, we would like to point out that our model can also be used to test a very 
interesting method for computing critical points that has be suggested recently This 
method is based on the assumption that the critical point of a system in d dimensions is 
located at the maximum of 


where T and are the transfer matrices for the d — 1 dimensional system and two coupled 
d — 1 dimensional system respectively. It has been shown that the above method is valid for 
self-dual models. However, for non self-dual models it remains a conjecture. Unfortunately, 
we found from our simulations that the critical point of the two layered Ising model is 
Pc = 0.2759(7) so that it does not appear to agree with the critical point obtained in Ref |Q. 

4 Critical Indices 


To establish the universality class, we need to investigate the hnite size scaling (FSS) ii 
behavior and determine the critical indices. 


For = ^2 = I//2 = /?, we determine u using 


dUL 

dp 


(1 - Ul) 


(<E> -2 


< mpE > 
< m? > 


< >\ 

> ) 


( 12 ) 


where E is the Ising energy per site. At Pc, dU^/dp should scale as As shown in 

Figure if we take Pc = 0.2759, we hnd u = 0.99(3). This is to be compared with the exact 
value V = 1 for the 2d Ising model. Next we measured the magnetization at the critical 
point Pc = 0.2759. Fitting to the scaling form 

m ~ , (13) 


as shown in Figure we get Pjv = 0.124(2). This should be compared to the exact value 
of 0.125 To get the anomalous dimension, rj, we measured the value of the susceptibility at 
Pc = 0.2759 as a function of L and htted to the form 

Xma. ~ . (14) 


As shown in Figure]^, we found r] = 0.27(6) in comparison to the exact value of 0.25. All the 
results on the critical exponents are summarized in Table The rest of the critical indices 
can be computed through scaling laws 0 a = 2 — i'd,'j/h' = 2 — ri,6 = {d + 2 — ri)/{d — 2+r]). 
Thus we conclude that all the critical indices of our model are the same as those of the 2d 
Ising model. Consequently these two models are in the same Universality class and they 
have to the same continuum limit: a 2d free Fermion system |I0|. 
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Figure 6: From the finite size scaling of X L'^-r one obtains the anomalous dimension 
ri = 0.27(6). 



Table 1: The critical exponents of the 2d Ising model and the coupled 2d Ising model. They 
seem to be identical with in the accuracy of our computation. 





5 Mass Spectrum 


We have established that the coupled Ising model in Eq. (|I]) and the 2d Ising model obey 
the same hnite size scaling laws and therefore belong to the same universality class. However 
the mass spectrum still remains a puzzle. Although there are two fields ax and Sx, we would 
like to know if there is one or two particles in the low energy spectrum. 


Let us introduce the Fourier transformation of the helds as 


<v)- ^d/2^^xe 

(15) 


(16) 

The propagators in Fourier space can be written as 


r(r.\ - ( < ^(P)^(-P) > < ^(p)s(-p) > \ 

^ y < s(p)cr(—p) > < s{p)s{—p) > j 

(17) 

Near p^ = 0, the two-point function can be expanded as 


G~^{p) = A + Bp^ + 0(p'^) , 

(18) 


where A and B are two-by-two matrices. B leads to a wavefunction renormalization. Thus 
the matrix M = B~^A is the mass matrix. The eigenvalues of M are the masses of the 
model (Here for simplicity, we dehne the masses at the zero momentum instead of the pole 
of the propagator). Numerically, we have found that it is fairly easy to fit G~}{p) 

near = 0 to a linear form, and G~l{p) is essentially independent of p^. We have obtained 
matrices A and B in Eq. (|T^ for y = 0.50 , K 2 = 0.20 and for various values of ki so that we 
cross the phase boundary. The results for mass eigenvalues are shown in Figure |^. It is very 
clear that while the small mass eigenvalue goes to zero at the critical point, the heavy mass 
stays at 0(1) in lattice unit. We also htted the light mass eigenvalue to the scaling form 

m? ~ (k — , (19) 

and found that for z/ = 1 one gets a very good fit of the data, (see Figure ||.) The deviation 
form the linear fit close to the critical point is due to FSS effects. This is another indication 
that the critical exponent v is indeed 1 as for the 2d Ising model. 

From the analysis of the propagator, we can learn two important things: (1) only one 
particle remains in the physical spectrum even though there appears to be extra degrees of 
freedom at the cut-off scale; (2) the physical held is a linear combination of bare helds ax 
and Sx- 


6 Conclusions 


In this paper we have studied a 2d spin model composed of two Ising systems interacting 
through a Yukawa coupling. There were two main objectives of our ehort. One was to show 
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Figure 7: The mass matrix eigenvalues at K 2 = 0.20 and y = 0.50. One goes to zero while 
the other remains finite at the eritieal point = 0.325. The mass unit is 1/a. 



K2 - 0.325 


Figure 8: The linear fitting of the low mass around the critical point confirms the fact that 
u is 1. The mass unit is 1/a. 
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that this form of coupling leads to the same universality class as a single 2d Ising model; 
the other was to show that only a single set of spin waves survives in the continuum limit. 
Our study of the critical exponents does support universality. All critical exponents agree 
to the accuracy of our computation. In addition the analysis of the two point correlation 
function supports the assertion that only one mass that goes to zero at criticality or that 
the continuum theory has a single mass scale. It is encouraging that this can be observed 
numerically in a two channel correlations function, when in reality the scalar excitations do 
not lie in a two dimensional subspace, but must include the multibody states. Nonetheless 
the freezing out of extra degrees of freedom is observable. Thus we were able to get solid 
numerical evidence that the continuum theory for the two spin Yukawa model is identical to 
the continuum limit of a single 2d Ising model. 

This result is not surprising from the viewpoint of renormalization group. One can 
think the integration over one set of spins as a “blocking” transformation which will induce 
irrelevant operators for the action of the 2d Ising model. As we noted in the introduction 
this result is almost “obvious”, but we are reassured by numerical “proof” and encourage 
to extend these numerical method to subtler coupled Higgs systems where our plausibility 
arguments might not be convincing. Our model does serves as an illustration for the key 
feature of Chirally Extended QCD Q] which, although involves extra degrees of freedom 
on the lattice, is expected to lead to the same continuum theory as the standard Wilson 
Lattice QCD. The numerical success of this demonstration gives us hope that more complex 
coupled systems can be studied by similar techniques. In the near future, we will attempt 
to extend these methods to 2d and 4d Yukawa coupled chiral spin models. 
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